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ABSTRACT 

We report on 8.4 GHz Very Long Baseline Interferometry (VLBI) observations of the type II-P supernova SN2004et in the spi- 
ral galaxy NGC6946, made on 20 February 2005 (151 days after explosion). The Very Large Array (VLA) flux density was 
1.23±0.07 mJy, corresponding to an isotropic luminosity at 8.4 GHz of (4.45±0.3)xl0-^ ergs"' Hz"' and a brightness temperature 
of (1.3±0.3)xlO** K. We also provide an improved source position, accurate to about 0.5 mas in each coordinate. The VLBI im- 
age shows a clear asymmetry. From model fitting of the size of the radio emission, we estimate a minimum expansion velocity of 
15,700±2,000kms"'. This velocity is more than twice the expected mean expansion velocity estimated from a synchrotron self- 
absorbed emission model, thus suggesting that synchrotron self-absorption is not relevant for this supernova. With the benefit of an 
optical spectrum obtained 12 days after explosion, we favor an emission model which consists of two hot spots on an underlying 
expanding shell of width comparable to that of SN 1993 J. 

Key words, galaxies: individual: NGC 6946 - radio continuum: stars - supernovae: individual: SN2004et 
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' 1. Introduction Li et al. ([2005') found a progenitor candidate for this super- 

C nova by analyzing Canada-France-Hawaii Telescope (CFHT) 

O Supernova SN2004et in galaxy NGC 6946 was discovered by images of the field made before the supernova explosion. The 

> ! Zwitter et al. (i2004|) on 27 September 2004. Yamaoka et al. progenitor seems to have been a yellow supergiant that may 

' (|2004|l reported data from the robotic Telescope a Action Rapide have previously experienced a red supergiant phase, but it could 

pour les Objects Transitoires (TAROT) telescope from which the also have consisted of a binary system of red and blue super- 

U supernova explosion date could be well constrained. Nothing giants. This latter case woul d impl y a progenitor similar to that 

brighter than magnitude 19.4+1.2 was seen on September for SN 1993J (Van Dvk et al. l20^ . 

22.0 1 7 and the supernova was detected at magnitude 15.17+0.16 

on September 22.983. Following Yamaoka et al., we assume Spectroscopic observations of Zwitter et al. (120041 1 and 

22.0 September 2004 as the explosion date. Filippenko et al. (2004) suggest that SN 2004et is a Type II su- 

NGC6946 is a starburst spiral galaxy with a high super- "ence, according to the emission model of Chevalier 

nova production rate. In fact, SN 2004et is the eighth supernova U9^, tested with other supernovae like SN 1993J (e. g., 

event registered in this galaxy since 1917. The distance esti- Marcaide et al. [1993, we expected that this supernova would 

mates to NGC 6946 range from 5.1 Mpc (de Vaucouleurs fT979] l ^mit detectable radio synchrotron radiation sometime after its 

up to 10.1 Mpc (Sandage & Tammann [TMB . Here we adopt oP^^al detection. Indeed, the radio emission from SN2004et 

5.5+0. 1 Mpc as the distance to NGC 6946 following Pierce was detected and momtored at 5 GHz with the Multi-Element 

( fT994l l. who estimated it using the Tully-Fisher HI relation. Most ^'"^^'^ Interferometer Network (MERLI N) and the VLA 

of the distance values quoted in the literature vary between 5.2 arrays from October 2004 onwards (Argo et al. 5005, Kelley et 

and 5.5 Mpc. Sahu et al. (2006J derived a distance of 5.6 Mpc ^1. m preparation). Momtoring was also carried out at 6.0 GHz 

to NGC 6946 as an average of the literature values, including the ^nd 1 .6 GHz by these authors. Unfortunately, the monitoring ob- 

result of their analysis for SN 2004et itself, which is in very good servations by Argo et al. were unable to constrain the date of 

agreement with the adopted 5.5 Mpc distance. *e 5 GHz peak emission, because scheduling consti-aints meant 

that no observations could be made between 3-17 November 

2004 (see Fig. 2 of Argo et al. 2005). From their observations. 

Send offprint requests to: I. Marti- Vidal Argo et al. ([2005) estimated the 5 GHz peak flux density as 

* Deceased ~2.5mJy in early November 2004. Kelley et al. (in prepara- 
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tion), using observations of SN 2004et at 6 different frequencies 
with the VLA, fitted the flux density evolution of the supernova 
to the "mini-shell" model of Chevalier (1982a 1982b) and ob- 
tained an estimate of the SN 2004et peak flux density at 5 GHz 
of 2.45+0.08 mJy on day 66+6 after explosion (assuming 22.0 
September 2004 as the explosion date). 

We arranged very-long-baseline-interferometry (VLBI) ob- 
servations of SN2004et in order to learn as much as possible 
about the radio emitting structure, the expansion scenario, and 
the emission mechanisms. In Sect. 2 we describe our observa- 
tions; in Sect. 3 we present our VLBI image of SN 2004et and 
the results of fits of several emission models to the data; in Sect. 
4 we discuss the magnetic field strength and the amount of en- 
ergy involved in the radio emission, together with the expansion 
velocity of the emitting region and the importance of possible 
absorption mechanisms; in Sect. 5 we present our conclusions. 



2. Observations and data reduction 

2.1. The VLBI data 

In our observations, made on 20 February 2005, the partici- 
pating stations were the complete Very Long Baseline Array 
(VLBA) (10 identical antennae of 25 m diameter spread over 
the US from the Virgin Islands to Hawaii), the phased VLA 
(equivalent to 130 m diameter. New Mexico, USA), the Green 
Bank Telescope (100 m, West Virginia, USA), Goldstone (70 m, 
Cahfornia, USA), Robledo (70 m, Madrid, Spain), Effelsberg 
(100 m, Bonn, Germany), and Medicina (32 m, Bologna, Italy). 

The recording rate was 256 Mbps, with 2-bit sampling and 
single polarization mode, covering a total synthesized band- 
width of 64 MHz (except at the phased VLA, where the ef- 
fective bandwidth was 50 MHz). These data were correlated at 
the Array Operations Center of the National Radio Astronomy 
Observatory (NRAO) in Socorro (New Mexico, USA) with an 
averaging time of 2 seconds. 

The 12-hour long, 8.4 GHz VLBI phase-reference observa- 
tions were designed to maximize the probability for the detection 
of supernova SN2004et. We used the source J2022+614 as the 
phase calibrator (a 3.0 Jy quasar at an angular distance of 2.2° 
from the supernova). We spent 25% of the observing time on the 
calibrator, with duty cycles of approximately 4 minutes. Every 
22 minutes we also observed a secondary phase calibrator, the 
quasar J2035+582, with a 0. 18 Jy total flux density and an angu- 
lar separation of 1.76° from the supernova. 

The correlation coefficients were then read into the NRAO 
Astronomical Image Processing System (aips) for calibration 
purposes. A manual phase alignment between the 8 different IPs 
was obtained by fringe fitting a selected scan of the calibrator 
source J2022+614. 

We then performed the amplitude calibration using gain 
curves for each antenna (information provided by NRAO and 
the European VLBI Network, EVN) and system temperatures 
measured by each station during the observations. 

Such a scheme is appropriate for all stations except the 
phased VLA. In that case, the amplitude calibration requires pre- 
vious knowledge of the flux densities of all sources. For consis- 
tency, we first mapped the calibrators using the VLBI data from 
all stations but the VLA. Then, we used the total fluxes obtained 
in the mapping of J2022+614 and J2035+582 to calibrate the 
phased VLA. 

The J2022+614 data were then fringe-fitted in a standard 
manner and the resulting calibrated visibilities were exported to 
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Fig. 1. Contour image of our phase calibrator source, 
J2022+614, at 8.4GHz. The dynamic range (peak/rms) obtained 
for this image is 68,000:1. The FWHM of the CLEAN beam is 
shown at bottom left. 



the Caltech program difmap (Shepherd et al. 11995 1) for imag- 
ing purposes using standard phase (and later also amplitude) 
self-calibration techniques. The phase corrections derived from 
self-calibration differed from the corrections derived from fringe 
fitting by typically 5- 10 degrees, and the amplitudes after self- 
calibration differed from the system temperature calibrated am- 
plitudes by typically 2-5%. The resulting map of J2022+614 
is shown in Fig.[T] Afterwards, the J2022+614 model obtained 
in DIFMAP was read back into aips for a second fringe-fitting in 
which we obtained residual phases, delays, and rate solutions 
nearly free of structure contributions. 

These fringe-fitted solutions were then linearly interpolated 
to the times of the SN 2004et scans and were used for the calibra- 
tion of the phases, delays and delay rates of the VLBI SN 2004et 
data. Finally, we exported the calibrated supernova data for hy- 
brid mapping in difmap. 

2.2. The VLA data 

The VLA provided two kinds of data: VLBI data from the 
phased array and individual antenna data. The correlation co- 
efficients for all the pairs of antennas could thus be obtained and 
used for wide-field imaging around SN 2004et and for obtaining 
an accurate estimate of the total flux density of SN 2004et. 

We also used J2022+614 as the flux calibrator for the VLA 
data. The radio image contours thus obtained are shown in Fig.|2] 
superimposed on an optical image. The map is affected by band- 
width smearing effects. The radio sources identified in Fig.|2]had 
previously been listed by Hyman et al. (2000). Since the correla- 
tion was centered on SN 2004et, its flux density was accurately 
determined as 1.23+0.07 mJy. 
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SN2004ET 8.4 GHZ WITH 1 CLEAN COMPONENT REMOVED 



Fig. 2. Radio image contours from the VLA superimposed on an 
optical image of NGC 6946 (taken with the TROBAR telescope 
of the Observatori Astronomic de la Universitat de Valencia). 
The FWHM of the CLEAN beam is shown at bottom left. 
Contours are at 0.02, 0.08, 0.32, 0.64, and 1 times the peak flux 
density, which is 1 .23 mJy/Beam. 
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Relative RA (mas) 
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3. VLBI image of SN2004et 

In DiEMAP, we deleted the bad visibility amplitudes of SN 2004et 
(i. e., less than 0.1% of the data) and Fourier inverted the re- 
maining data into the sky plane. We then applied the natural 
weighting scheme to the visibilities for sensitivity optimization 
and detected a nearly point emitter (showing a small elongation 
in the South-East direction) whose peak was displaced 5.1 mas 
West and 14.7 mas South from the coordinates at which we 
had centered the correlation (we used the coordinates given by 
Beswick et al. 2004). The emitter is undoubtedly SN2004et. 
Thus, a better position of the supernova, based on a phase- 
reference to J2022-1-614, a source of the International Celestial 
Reference Frame (ICRF), is: a = 20^ 35™ 25.3596^ and 5 = -1-60° 
07' 17.7203" (equinox J2000.0, ±0.5 ma^] in each coordinate). 

Given the very low flux-density of SN2004et, we did not 
apply any phase self-calibration to the data. We fitted several 
emission models to the real and imaginary parts of the visibilities 
using the difmap task modelfit. The errors of the fitted parameters 
of the models correspond to the statistical uncertainties taken 
from the covariance matrices of the fits. A point source model 
fitted to the data has a total flux density of 0.63+0.03 mJy. After 
subtraction of this point source model from the data, there is still 
clear extra emission at the level of 0.24 mJy beam"', indicating 
structure in the supernova emission. A map of this extra emission 
is shown in Fig.|3] 

This extra emission could be due to a bad antenna calibra- 
tion, so we checked on this possibility by making a map of the 
source leaving out one antenna of the array at a time. We also im- 
aged the first and second halves of these data separately to check 
the structure repeatability. All the resulting maps showed the ex- 
tra emission when we subtracted a point source from the uv-data. 
Thus, we conclude that this emission belongs to SN 2004et and 
is not produced by imperfect calibration. 

Using two point emitters to fit the data, we recovered a 
flux density of (0.57+0.05) mJy for the first component and 
(0.26+0.05) mJy for the second, located at an angular distance 



' This error corresponds to the beamwidth of the VLBI interferome- 
ter applying the natural weighting scheme to the SN2004et visibilities. 



Fig. 3. Residuals after the subtraction of a single point source to 
SN 2004et (see text in Sect. 3). 



of 0.5+0.06 mas and at a position angle of (136+4)° with respect 
to the strongest component. 

Given that the total recovered flux density in the uv data us- 
ing CLEAN adds up to 0.87 mJy (5% higher than the flux den- 
sity fitted), and the VLA flux density estimate is 1 .23 mJy, there 
must be some missing flux density in our VLBI data. Most of 
this missing flux density could be due to a loss of coherence 
in the phase-referenced VLBI data compared to the VLA data. 
This part of the missing flux could be recovered using phase 
self-calibration, were the flux density of the source high enough. 
However, given that this is not the case, we decided to leave the 
phases unchanged and consider this lack of flux density as an 
instrumental limitation. 

Disks or shell-like sources could also be used to fit the data. 
Indeed, we show the results of some of these fits in Table [1] We 
find that for the fit to two point emitters the resulting rms of the 
residuals is slightly (but not significantly) smaller than for any 
of the other fits, which are very similar. The similarity between 
the rms of the fits of all the diff'erent emission models indicates 
that the source is only barely resolved by the interferometer and 
it is not possible to distinguish between more detailed models. 

We show the VLBI image of SN 2004et in Fig. |4] The peak 
brightness is 0.64 mJy beam ' and the rms of the residuals 
0.03 mJy beam"'. 



4. Discussion 

4.1. Supernova energy and magnetic field 

The SN 2004et flux density measured by the VLA on February 
2005 using J2022+614 as the flux density calibrator is 
1 .23 ±0.07 mJy. For a distance to NGC 6946 of 5 .5 ±0. 1 Mpc, this 
flux density corresponds to an isotropic luminosity at 8.4 GHz 
of (4.5±0.3)xl0^^ erg s"' Hz"'. The corresponding brightness 
temperature, if we use the size determined from the fit to two 
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Table 1. Results of different models fit to the SN 2004et data. 

' The assumed shell width is 30% of the outer radius. 

^ The diameter refers to the separation of the point emitters. 



Model 




Results 






Diameter 


Total flux density 


rms of residuals 




(mas) 


(mJy) 


(/ijy/beam) 


Uniform disk 


0.62 + 0.08 


0.81 ± 0.05 


27 


Optically thin sphere 


0.8 + 0.1 


0.83 ± 0.05 


27 


Optically thin shell' 


0.70 ± 0.08 


0.82 ± 0.05 


27 


Two point emitters^ 


0.50 ± 0.06 


0.83 ± 0.07 


26 



SN2004ET AT 8.4 GHZ 




4 3 2 1 -1 -2 -3 

Relative RA (mas) 
Center at RA = 20h 35m 25.35921749s , Dec = 60d 07' 17.7203001" 
Peak flux = 6.3860E-04 JY/BEAM 

Levs = 6.386E-05 * (-1 , 1 , 2, 3, 4, 5, 6, 7, 8, 9 1 0) JY/BEAM 



Fig. 4. Contour map of SN2004et obtained from our VLBI 
data by model fitting two point sources (see text in Sect. 3). 
The lowest positive level shown is 10% of the peak, which is 
0.64mJy/Beam. The FWHM of the CLEAN beam is shown at 
bottom left. 



point emitters, is (1.3+0.3)xlO^ K, indicating a non-thermal syn- 
chrotron origin of the radio emission. 

We did not obtain spectral information of the radio emission 
coming from SN2004et in our VLBI observations. However, 
Kelley et al. (in preparation) observed SN 2004et with the VLA 
at several frequencies (ranging from 1.4 GHz to 43.3 GHz) dur- 
ing several hundred of days following the supernova explosion. 
These authors fit the evolution of the resulting radio light curves 
to a given parametric model, from which they estimate a spectral 
index, a = -0.87 + 0.04, for the radio emission in the optically 
thin phase of the supernova evolution. Given that the supernova 
was in its optically thin phase during our VLBI observations, 
we used the spectral index estimated by Kelley et al. From the 
previously derived luminosity at 8.4 GHz, we estimate a total 
integrated radio luminosity of (1.3+0.1)xlO-'^erg s"' between 
1.4 GHz and 43.3 GHz. 

Using this result together with the SN 2004et size estimate 
from our VLBI observations, we can derive the minimum total 
energy in the radiating region of the supernova and the corre- 



sponding equipartition magnetic field. From Pacholczyk ( 119701 ) 
(see chapter 7, formulae 14 and 15) we have; 

£n,i„=CB(l+fc)'*^V^^V/^Lr (1) 

Beq = (4.5ci2(l + k)/4>f'' R-^'^Lf (2) 

where ci2 and cn are functions of the spectral index, a, and the 
minimum and maximum frequencies considered in the spectrum 
integration, (f> is the filling factor of the emitting region, R is the 
source radius, Lr is the integrated radio luminosity, and k is the 
ratio between the heavy particle energy and the electron energy. 
The mechanism of generation is not well known, hence k can 
vary from 1 to nip/me ^ 2x 10^. Thus, these formulae give us a 
range of possible values of the equipartition energy and magnetic 
field instead of an estimate. Also we do not know the value of the 
filling factor, cp (the fact that the emission may not be completely 
isotropic will also affect the estimate), but the effect of <p on the 
value of Emin and Bgq goes aproximately as the square root and 
the cube root, respectively, as can be seen in equations [1] and |2] 

Assuming that the radio emitting region corresponds to a 
shell of width equal to 30% of the outer radius (e. g., Marcaide 
et al. [1995 2007 ), we find that the filling factor, 4>, is 0.66, and 
that the minimum total energy (particles plus fields) of the emit- 
ting region (with a radius of 0.70+0.08 mas) ranges between 
8.2x10''^ erg and 4.22x10"*^ erg. The upper part of the range still 
represents less than one part in a thousand of the supernova emis- 
sion energy, ~10^' erg. 

The corresponding equipartition magnetic field ranges be- 
tween 37 mG and 260 mG. This range is similar to those ob- 
tained for other type II supernovae, like SN 1986 J (Perez-Torres 
et al. l"2002l l. SN 1993J (e. g., Perez-Torre s et al 1200 II I or, more 
recently, SN2001gd (Perez-ToiTes et al. |2005] l. Values in this 
range are very high compared to the typical (about 1 mG) mag- 
netic fields in the circumstellar medium of a massive star Given 
that the hydrodynamic compression from the supernova shock 
only increases the circumstellar magnetic field by a factor of 
about 4 (e.g., Dyson & Williams ,1980) , other field amplifica- 
tion mechanisms, i. e. turbulent amplification possibly due to 
Rayleigh-Taylor instabilities (Chevalier 1 1 982bl l, must be taking 
place in the emitting region of SN 2004et. 

4.2. Supernova structure 

Our results show that the emission structure of supernova 
SN2004et is not spherically symmetric. The residuals of the fit 
of a point emitter to the VLBI data, shown in Fig.|3] indicate that 
the emitting structure has an angular asymmetry. 

The structure of the ejecta of type II-P supernovae are sup- 
posed to be very spherical, given that their extended hydro- 
gen envelope should "spherize" the possible asymmetries in the 
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shock (e. g., Khokhlov, Hoflich & Wang |200TT l. Thus, the asym- 
metry found in SN 2004et is Ukely to come from the structure of 
the circumstellar material. 

Assuming a circumstellar origin for the asymmetry of the 
SN 2004et emission, we may interpret such asymmetry in two 
ways: 

Model 1. A shell with two hot spots, one on each "side" 
of the emission of the shell. These hot spots could be due to 
anisotropics of the circumstellar density distribution (probably 
due to asymmetries in the progenitor pre-supernova wind) and/or 
anisotropics of the magnetic field. This kind of anisotropic shell 
has been found in the case of SN 1986J (Baitel et al. I199U 
Perez-Torres et al. 2002, Bietenholz et al. .2004) . The fact that 
we recovered a much lower flux density with VLBI than with 
the VLA is consistent with this model. The two hot spots of 
the radio structure would emit the major part of the total flux 
density of the source, and the emission from the shell, being 
faint, would be buried in the thermal noise of the map. A flux- 
density distribution with a peak surface brightness just under 
0.1 mJy beam ' (peak flux density of the noise; see the spread 
contours of Fig. [3j distributed over an area of 2 VLBI beams 
(roughly, the area of the radio structure), would account for a 
flux density of ~0.2 mJy, that is, half of the difference between 
the flux density recovered by the VLA and VLBI. 

Model 2. An expanding shell plus a protrusion. This pro- 
trusion could be produced by an asymmetry in the circumstel- 
lar medium distribution, like an axi-symmetric region with a 
lower circumstellar density. Such asymmetry could result in a 
faster expansion of the portion of the forward shock located in 
that region. According to hydrodynamical simulations (Blondin, 
Lundqvist, & Chevalier 1996), the length of the protrusion could 
be more than twice the radius of the main shell, but the time re- 
quired for the formation of such a structure in the forward shock 
would be of the order of several years, much larger than the age 
of SN 2004et at the time of our observations. 

Given the youth of the supernova (151 days at the epoch of 
our observations), the size of the protrusion should be much 
smaller than the mean size of the shell. This argument does not 
favour Model 2 as the origin of the observed structure. 

4.3. Expansion velocity 

The different models shown in Table [1] yield different estimates 
of the size of SN 2004et. We do not know which one of these 
models (or another) best represents the true emitting structure 
of SN2004et. Thus, we will estimate the minimum expansion 
velocity compatible with our data using the smallest of our size 
estimates. This separation (see Table [1) is 0.50+0.06 mas and 
corresponds to the model of two compact emitters. 

Assuming a distance of 5.5+0. 1 Mpc for NGC 6946 and 22.0 
September 2004 as the explosion date of the supernova, we find 
that the radial distance between the point components translates 
into a mean separation speed of 31, 400+4,000 km s"' between 
them. Assuming also that the two compact emitters come from 
two diametrically opposed hot spots on an expanding shell, the 
center of the explosion would fall half way between them and the 
expansion velocity of the shell would be 15,700±2,000kms"'. 
This velocity is high compared with the theoretical predictions 
of a synchrotron self-absorbed model, as will be shown below. 

If we use the size estimates of the other models in Table [T] 
the derived expansion velocities are higher than 15,700kms 
Fitting the data to a disk with uniform emission (unrealis- 
tic, because the circumstellar region is not optically thick 
at this epoch), yields an expansion speed of the ejecta of 
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Fig. 5. Red part of an optical spectrum of SN 2004et obtained 
at the David Dunlap Observatory (see Vinko et al., in prepara- 
tion). The spectrum is normalized to the continuum. The vertical 
dashed line marks the position of the blue absorption minimum 
of Ha, corresponding to the velocity of the Ha "photosphere" 
(see text). 



19, 500+2, 500 km s ' ; fitting to an optically thin shell-like model 
yields a speed 22, 000+2, 500km s"', and fitting to an optically 
thin filled sphere results in a speed of 25,200+3,000 km s"' . 

An early optical spectrum of SN 2004et is shown in Fig. |5] 
and shows characteristic P Cygni profiles of an expanding en- 
velope. This spectrum was taken on 3 October 2004 (7 days af- 
ter discovery and 12 days after explosion) at the David Dunlap 
Observatory in Canada with the Cassegrain spectrograph on the 
74 inch (188 cm) telescope as part of observations that will be 
fully described by Vinko et al. (in preparation). From the posi- 
tion of the blue absorption minimum of the Ha line, an expan- 
sion velocity of the line of 11,900 + 600km s"' is estimated, 
once the redshift of the host galaxy is corrected. This velocity 
represents the velocity of the Ha "photosphere", i. e. where the 
optical depth becomes r ~ 1. Because at early phases the ejecta 
are almost fully ionized, this layer should be close to the top of 
the ejecta. Thus, the velocity from the Ha line should also be 
close to (although somewhat lower than) the expansion velocity 
of the ejecta. 

Sahu et al. ( 120061 1 have reported expansion speeds from op- 
tical observations starting on day 30 after explosion (25 days 
after discovery). They identify a low velocity component in Ha, 
which is the same as that shown in Fig.|5]but for later days and 
hence at lower velocities. They also identify a high velocity com- 
ponent for days 30 through 90. Such high velocity feature is ex- 
pected to arise from the interaction between the CSM and the 
SN shock wave (Chugai et al. l2007l l. It corresponds to a feature 
placed on the absorption side of the P Cygni profile of the lower 
velocity Ha, both in their spectra and in the spectrum shown in 
Fig E] However, this component appears very weak and almost 
imperceptible in Fig. |5] We are not convinced that it provides a 
good estimate of the early expansion speed and we will not use 
hereafter Also, Chugai et al. (I2007I I did not find convincing ev- 
idence for the high-velocity Ha feature in the middle and late 
plateau phases, when such features should be more pronounced. 

According to the model of Chevalier ( |1982bl ), the opti- 
cal emission from young supernovae originates in the shocked 
ejecta region, located behind the shocked circumstellar region. 
This implies that the maximum velocities inferred from the op- 
tical lines should always be smaller than the velocities inferred 
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from the VLBI observations of the radio emitting structure. The 
latter radio emission would originate in the shocked circumstel- 
lar region. In fact, according to the self-similar expansion model 
of Chevalier (|1982al [1982b ), the ratio between the maximum 
velocity inferred from the optical spectral lines in the very early 
days and the expansion velocity of the radio shell should be sim- 
ilar to the ratio between the inner and the outer radii of the shell- 
like radio emitting structure. Marcaide et al. (I2007I I find a ratio of 
0.69+0.02 between the inner and the outer radii of the shell-like 
structure of supernova SN 1993J. 

For SN 2004et we find that the ratio between the maximum 
velocity inferred from the optical spectral lines in the very early 
days and the expansion velocity of the radio shell is about 0.76 
for the model yielding the smallest expansion speeds from the 
radio data. Thus, if the components of the two-point model are 
hot spots on the outer part of an underlying shell (which we do 
not detect, as explained above, due to its faintness) the width of 
such a shell would be comparable to that of SN 1993J. In aiTiv- 
ing at this conclusion, we have not considered the deceleration 
of the supernova ejecta between the epoch of the optical obser- 
vations and the epoch of the VLBI observations. This neglect 
implies that the optical line velocity at our VLBI epoch could 
be lower than 11, 900 km s"' . Thus, also taking into account that 
15,700km s"' is the smallest expansion speed estimated from 
these radio data, we conclude that the real ratio between the 
inner and the outer radii of the supernova radio shell could be 
somewhat smaller than 0.76. 

For cases of supernovae where synchrotron self-absorption 
(SS A) dominates over the circumstellar free-free absorption, and 
where there is equipartition of the energy between relativistic 
particles and magnetic fields. Chevalier ( 1998 ) found a relation- 
ship between the observing frequency, the size of the emitting 
region, the peak intensity of the supernova light curve, the su- 
pernova age at the peak, and the spectral index of the energy 
distribution of the relativistic electrons. For such a case, one can 
estimate the expansion velocity of a supernova based on mea- 
surements of the peak luminosity at a given frequency, the su- 
pernova age at the peak, and the synchrotron spectral index, as 
follows; 

/ p \0-4729 

yssA(kms-i) ^ 5 37gg ^ io6(^0)-''"54i I _Ll x 



IMpc/ Vl.OGHz/ \lday/ 

where fi is the ratio of the relativistic particle energy density 
to the magnetic field energy density (since we assume energy 
equipartition, jS - -^^y^, Pacholcz vk 1 1 9701 chapter 7), (f> is the 
filling factor of the emitting region for a sphere (we assume a 
shell of width equal to 30% of the outer radius, which yields 
(p - 0.66), F-p is the flux density at the peak, D is the distance, v is 
the observing frequency, and tp is the age at the peak. We use the 
synchrotron spectral index found by Kelley et al. (in preparation) 
to set the index of the electron energy distribution to y = 2.74 
(the energy distribution of the relativistic electrons is assumed to 
be of the form oc and y - \ - 2a). 

When we use equation [3] and assume for Fp the 5 GHz peak 
flux density published by Argo et al. ( I2005I I (i. e., Fp - 2.5 mJy 
at tp - 42 days), we estimate a SSA expansion velocity of the 
radio shell of SN2004et of 9,700+l,800kms-'. This is much 
lower than the estimate for each of the models in Table [T] If 
instead we use the 5 GHz peak flux density fitted by Kelley et al. 



(in preparation) at day 66+6 after explosion, which appears more 
accurate given the extended time coverage of their observations 
near the emission peak at 5 GHz, the SSA expansion velocity 
of the radio shell is only 6,100±l,100kms"'. This result is less 
than half of any of the estimates based on our VLBI data. 

As noted by Chevalier ( 119981 ), if the free-free absorption 
from the circumstellar medium is comparable to, or even higher 
than, the synchrotron self-absorption, the peak of emission will 
be defined by the free-free absorption evolution. This will pro- 
duce a negative bias in the estimate of the ejecta expansion ve- 
locity, as we seem to find for SN 2004et. Comparison of VLBI 
results and radio light curves for other supernovae also gives 
important negative biases (see Chevalier ,1998 ), but at present 
there is an insufficient number of cases (in fact, only SN 1979C, 
SN 1986 J, and SN 1993 J) to reach firm conclusions about the 
typical values of these biases. For supernovae SN 1986J and 
SN 1979C, the bias ratios Vvlbi/ Vssa (-2.56 and -2.60, respec- 
tively) are similar to that for SN 2004et (if we use the minimum 
size compatible with our data), but for SN 1993 J this bias ratio 
is smaller (~1.5). 

4.4. Comparison with other supernovae in NGC 6946 

From those eight supernovae known to have exploded in 
NGC 6946, only four have been detected as radio supernovae: 
SN1968D, SN1980K, SN2002hh, and SN2004et. For two 
(SN 1980K and SN20004et), it has been possible to fit the light 
curves (Weiler et al. 119921 Kelley et al. in preparation) in terms 
of the "mini-shell" model (Chevalier 1982al[T982b] l. 

There are a number of similarities between SN 1980K and 
SN2004et. The flux-density power-law time decay parameter 
takes the values -0.62 and -0.74, respectively, and the power- 
law parameter of the evolution of the external opacity with time 
takes the values -2.46 and -3.07, respectively. But there are also 
significant differences between these supernovae: the spectral in- 
dex is -0.52 for SN 1980K, while for SN 2004et it is -0.87; the 
time after explosion needed to reach the peak spectral luminos- 
ity was 134+2 days for SN 1980K and is only 66+6 days for 
SN 2004et. For the case of SN 1980K (Weiler et al. 1992 ), inter- 
nal and clumpy opacity contributions are not required to obtain 
a good fit, and the light curves are well characterized by an opti- 
cally thin synchrotron radio emission absorbed by an uniform 
external mediurrQ. However, for SN2004et (see Kelley et al. 
in preparation), the absorption produced by a uniform external 
medium cannot fit the data, but a good fit is obtained for absorp- 
tion by a clumpy or filamentary external medium. The presence 
of this clumpy external medium could explain the asymmetry in 
the radio structure that we have found in our VLBI image. On 
the one hand, an expansion into a clumpy medium would yield 
an irregular interaction between the ejecta and the circumstel- 
lar medium and, thus, the evolution of the expanding material 
would be asymmetric. On this account. Model 2 would appear 
supported. On the other hand, the synchrotron emission, being 
proportional to the relativistic electron density in the circumstel- 
lar region, would be larger in the clumps. Thus, these clumps in 
the CSM would give rise to several hot spots in the expanding 
radio shell and would appear to give support to Model 1. The 
asymmetry found in our VLBI image is, then, consistent with 
the expansion of the supernova ejecta into a clumpy medium. 



' For SN 1968D, it has not been possible to fit the light curves due to 
the late detection (at an age of ~26 years), but the measured flux densi- 
ties are consistent with the extrapolated radio behaviour of SN 1980K. 
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5. Conclusions 

We report on 8.4 GHz VLBI observations of supernova 
SN2004et in NGC 6946. The observations, made 151 days af- 
ter explosion, are the first VLBI observations of a type II-P su- 
pernova, and successfully detect the radio emission and yield an 
image of the source. This image shows evidence of structure. 

This structure has a geometry that could be interpreted either 
as emission coming from a shell with two hot spots {Model 1) or 
emission coming from a region expanding anisotropically and 
developing a protrusion (Model 2). However, such a protrusion 
would be too large for so young a supernova, to be compatible 
with the hydrodynamic simulations of Blondin, Lundqvist, & 
Chevalier (1996). 

The asymmetry found in our VLBI image can be explained 
if the supernova expands in the clumpy circumstellar medium 
suggested by the fits to the radio light curves of this supernova 
(Kelley et al. in preparation). 

The minimum expansion speed that we obtain is 
15,700±2,000kms"' for a two point emitter model. Using 
other models yields larger expansion speeds. Maximum velocity 
estimates obtained from optical spectral line observations give 
lower values than the velocities inferred from our VLBI data 
for all cases. Such differences between maximum optical and 
radio expansion speeds can be understood in the framework of 
the self-similar expansion model of Chevalier ( |1982al |1982bl ). 
In particular, the ratio of the maximum early Ha line velocity 
(corresponding in the Chevalier model to the maximum velocity 
of the ejecta) to the expansion velocity estimated from VLBI for 
a two point model is ~0.76, similar to that for SN 1993J. This 
ratio suggests that Model 7 is a plausible model. 

In any case, the expansion speed derived from our VLBI data 
is more than twice larger than the estimate obtained using the 
synchrotron self-absorbed model proposed by Chevalier (I1998I I. 
This discrepancy indicates that the synchrotron self-absorption 
must be negligible (i. e., small or nonexistent) compared to the 
circumstellar free-free absorption in SN 2004et at our observing 
epoch. 

The flux density obtained from our VL A data (1.23 mJy) and 
our minimum VLBI size estimate (diameter of 0.5 mas) trans- 
late into a brightness temperature of (1.3+0.3)xlO'*K, indicat- 
ing a synchrotron origin for the radiation. Using the spectral in- 
dex a - -0.87 + 0.04, determined by Kelley et al. (in prepa- 
ration) with data between 1.4 and 43.3 GHz, we find that the 
amount of energy injected into the emitting region is in the range 
(8-420)xl0'*^ erg, a small fraction (lower than 0.001) of the typ- 
ical emission energies in supernova explosions. The equiparti- 
tion magnetic field falls within (37-260) mG, indicating the ex- 
istence of highly efficient field amplification due perhaps to tur- 
bulent instabilities in the emitting region. 

Further VLBI observations of this supernova are not possi- 
ble given its low flux density, now well below the VLBI sensitiv- 
ity limit at the possible data rates at present. VLBI observations 
of similar. Type II-P, supernovae will be essential for the un- 
derstanding of the details of the radio emission of this kind of 
explosion. 



Observatori Astronomic of the UV for obtaining the optical image of NGC 6946. 
We thank the anonymous referee for his/her useful comments. 
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